The ^"^Sr optical frequency standard at PTB 
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With ^'^Sr atoms confined in a one dimensional optical lattice, the frequency of the optical clock 
transition 5s^ ^So - 5s5p ^Po has been determined to be 429 228 004 229 872.9(5) Hz. The transi- 
tion frequency was measured with the help of a fs-frequency comb against one of PTE's H-masers 
whose frequency was measured simultaneously by the PTB Cs fountain clock CSFl. The Sr optical 
frequency standard contributes with a fractional uncertainty of 1.5 x 10~^® to the total uncertainty. 
The agreement of the measured transition frequency with previous measurements at other institutes 
supports the status of this transition as secondary representation of the second with the currently 
smallest uncertainty. 

PACS numbers: 06.20.fa, 06.30.Ft, 32.30.Jc, 37.10.Jk 



INTRODUCTION 

The development of frequency standards, both, in the 
microwave [l|43| and optical [^421 domain is a very active 
field of research. The investigations reach out far beyond 
the aspects of metrology into questions of fundamental 
research [8-12|. The development of optical clocks using 
optical transitions as reference has provided a number 
of systems l^, ]A\ that have a stability far better than 
today's Cs fountain clocks. Several ones have already 
the capability to achieve a systematic uncertainty better 
than the best realizations of the second [1, H, 0| ■ 

The field of optical clocks can be divided into frequency 
standards with single ions (or as proposed few ions) and 
with a large number (~10'*) of neutral atoms. Single ion 
standards lead the race for the highest accuracy Q . Stan- 
dards with neutral atoms, however, promise an extraordi- 
nary stability due to the larger number of absorbers that 
are simultaneously interrogated. Among these, systems 
using the fermion ^^Sr as reference are widely distributed 
and are investigated to our knowledge in at least eight 
laboratories. 

The newly developed optical clocks with high stabil- 
ity allow measuring frequency ratios without any refer- 
ence to the realization of the second by Cs clocks. With 
stabilized optical links lB"iZ| these measurements can 
in principle be performed between distant clocks. How- 
ever, they do not yet connect optical frequency standards 
based on the same transition in different institutions 
to test their accuracy beyond the limit set by primary 
clocks. In this respect, only a few high accuracy tests 
of such standards have been reported so far 0, 13| ; and 
these have been within one lab. Frequency measurements 
of optical frequency standards against primary clocks are 
therefore of high interest to evaluate the performance of 
distant systems. 

So far, for ^^Sr lattice clocks high accuracy mea- 
surements have been performed in three laboratories 
17 . [Til . [20l |. All results show an excellent agreement. 



which has resulted in the acceptance of the 5s^ ^Sq - 
5s5p '^Fq clock transition as secondary representation of 
the second with currently the smallest uncertainty 13]. 

We present here the first frequency measurement of 
this transition performed with the Sr optical frequency 
standard at the Physikahsch-Technische Bundesanstalt 
(PTB). An uncertainty of 0.5 Hz or in fractional units of 
1.0 X 10~^^ has been achieved, confirming the previously 
published results. 

The paper is organized as follows: In Sec. □ the ex- 
perimental setup and interrogation sequence is outlined 
in detail. Systematic shifts of the clock transition and 
their uncertainty are discussed in Sec. □ followed by the 
frequency measurement and its results in Sec. „ 



EXPERIMENTAL SETUP 

Part s of th e experimental setup have been described 
before 
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-23|. For completeness, we give a detailed de- 
scription of the frequency standard as used to perform the 
absolute frequency measurement: *^Sr atoms are evapo- 
rated from a furnace at a temperature of 520 °C, col- 
limated by a 2 mm diameter pinhole. They are laser 
cooled in a Zeeman slower and then deflected and col- 
limated by a two dimensional optical molasses operated 
with light red detuned by 385 MHz (slower) and 10 MHz 
(molasses) from the 461 nm transition ^Sq - ^Pi (for 
the relevant atomic levels see Fig. [T]). The deflection of 
the cooled atomic beam into the capture region of the 
magneto-optical trap (MOT) ensures that no direct line- 
of-sight exists between the hot furnace and the trapping 
region. Additionally, the room-temperature pinhole be- 
tween furnace and Zeeman slower reduces the indirect 
blackbody radiation flux in the direction of the MOT. 

The MOT with a magnetic field gradient of 7.4 mT/cm 
is operated with light detuned by 32 MHz below the 
461 nm transition (natural line width: 32 MHz). The 
1 /e^ beam radii of the MOT beams are 5 mm and they 
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FIG. 1: Partial level scheme of Sr with the two cooling 
transitions (461 nm and 689 nm), repumpers (679 nm and 
707 nm), and the clock transition at 698 nm. The hyper- fine 
structure splitting of the '^Pi level is indicated on the right. 



have peak intensities of about 3.5 mW/cm^ each. In the 
MOT a cloud of ®^Sr atoms is cooled within 300 ms to 
about 3 mK. All laser beams are derived from a frequency 
doubled diode laser system whose frequency is stabilized 
to a long-term stable reference resonator. To increase 
the atom number, repump laser beams resonant with the 
707 nm - ^Si and 679 nm ^Pq - ^Si transitions are 
applied. To cover the hyperfine manifold of the ■^P2 level 
the 707 nm repump laser is rapidly scanned over an in- 
terval of 5 GHz. Its average frequency is stabilized with 
a commercial wavemeter. The 679 nm repump laser is 
frequency stabilized to a reference resonator. 

The pre-cooled atoms are then reloaded into a MOT 
with lowered magnetic field gradient of 0.7 mT/cm op- 
erated on the 689 nm intercombination line ^Sq - ^Pi 
with a natural linewidth of 7 kHz. The laser is tuned 
1.6 MHz below the resonance of the AF = F' — F = +1 
transition (with the total angular momentum F) and is 
frequency modulated at 30 kHz with a peak-to-peak ex- 
cursion of 3 MHz to increase the velocity capture range. 
The MOT beams have 1/e^ beam radii of 2.6 mm and 
peak intensities of about 5.5 mW/cm^. Beams from a 
second laser (the so-called stirring laser) are superim- 
posed with the cooling laser beams. The beam intensity, 
detuning, and frequency modulation with respect to the 
689 nm AF — transition are similar to the one of the 
cooling laser. Stirring light is applied to redistribute the 
population in the Zeeman sublevels and thus enhance the 
cooling efficiency, which is otherwise reduced by optical 
pumping to untrapped levels. This is due to strongly dif- 
ferent g-factors of the ^Sq and '^Pi levels of the cooling 
transitions that cause large variations in the detuning of 
the cooling beams from the Zeeman transitions (and thus 
excitation probabilities) that are not compensated by the 
Clebsch- Gordon coefficients f24!|. After a 90 ms long cool- 



ing phase a temperature of about 50 /iK is reached. The 
temperature is further lowered within 50 ms to ~5 /iK 
applying only the cooling laser beams without frequency 
modulation at reduced detuning (400 kHz) and intensity 
(440 iiW/cw?). 

A horizontally oriented and linearly polarized ID op- 
tical lattice is overlapped with the atomic cloud during 
the whole cooling sequence. The lattice is formed by 
a retro-reflected laser beam from a Ti:Sapphire laser at 
the magic wavelength at 813 nm [25|. The frequency 



of the lattice laser is stabilized to a fs frequency comb 
and monitored with a wavemeter with 1 MHz accuracy 
to unambiguously detect mode hops of the lattice laser 
followed by relocking to a different fs-comb mode. The 
lattice waist radius is about 32 //m; the trap depth is 
typically Ut — 125 iJ^, where = h?k^/2m is the pho- 
ton recoil energy with the wavevector k = 27r/A of the 
lattice light and m is the atomic mass of ^^Sr. Motional 
sidebands in the lattice are observed with frequencies of 
i^z ~ 77 kHz in axial and Vr ~ 440 Hz in radial direction. 

After turning off the light of the MOTs the atoms 
in the lattice have a temperature similar to the one in 
the last MOT phase. The atoms in the lattice are then 
spin-polarized to either the rap ~ +9/2 or m^? = —9/2 
^So Zeeman level by an appropriately circularly polarized 
light pulse of 1 ms length resonant with the AF = ^Sq - 
■^Pi transition. To define the quantization axis, the MOT 
quadrupole coils are used in approximately Helmholtz 
configuration to generate a homogeneous magnetic field 
of about 23 /iT. The magnetic field axis is parallel to the 
lattice polarization. After spin-polarization, the mag- 
netic field is increased to 1.8 mT and an intense tt pulse 
of around 4.8 ms duration, resonant with the ^Sq - ^Pq 
mp = +9/2 or = —9/2, Amp = clock transition, 
is applied. The clock laser beam is aligned with the axis 
of the optical lattice with a residual angle of less then 
0.2 mrad. It has a 1/e^ radius of about 65 /im and has 
linear polarization that is parallel to the orientation of 
the magnetic field and the linear lattice polarization. 

The light is delivered via a polarization maintain- 
ing optical fibre from an extended cavity diode laser in 
Littmann configuration. The laser frequency is offset by 
an acousto optical modulator (AOM) and is stabilized to 
a resonator made of ultra low expansion glass [21|. The 
resonator has a finesse of about 330 000 and is housed 
in a gold plated passive heat shield in an ultra high vac- 
uum system that is actively temperature stabilized. The 
whole reference cavity setup is mounted on a passive vi- 
bration isolation platform. A laser linewidth of 1 Hz was 
observed by comparison with other laser systems via a 
fs-laser comb 
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We use an active fibre noise cancellation system to 
stabilize the optical path length between the clock laser 
and the retro-reflection mirror of the optical lattice on 
which the clock laser beam is overlapped with the lat- 
tice. This mirror determines the position of the atoms. 
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FIG. 2: High resolution spectrum of one Zeeman component 
of *^Sr excited with a 90 ms spectroscopy pulse. A Fourier 
limited linewidth of 9 Hz and an excitation probability of 
about 85% is achieved. Each point is a single measurement 
requiring 0.625 s. 

By stabilizing the path length we ensure a constant phase 
of the clock laser field at the position of the atoms. 
The atoms are interrogated by pulsing the light deliv- 
ered through the fibre with an acousto-optical modula- 
tor (AOM). Thus, the stabilization must be operated in 
a pulsed mode, too. Special care has been taken to avoid 
transient phase chirps during switching (SeCo details will 
be presented elsewhere [27]). 

Having prepared typically 2 x 10^ ®^Sr atoms in a de- 
fined mp level of the ■^Pq state by the tt pulse described 
above, atoms remaining in the ^Sq state are pushed away 
by radiation pressure from a resonant 461 nm laser beam. 
Then the Amp — transition starting from the prepared 
extreme mp level is interrogated in a homogeneous field 
of 23 iiT with a 7r-pulse of typically 90 ms length. We 
observe Fourier limited spectra with high (de-)excitation 
probability (Fig. To determine the transition prob- 
ability we, first, detect ^Sq atoms by exciting the ^Sq 
- ^Pi transition and detecting the fluorescence with a 
photo diode. Then the ground state atoms are removed 
and atoms in the '^Pq state are pumped via the '^Si state 
to the ^So ground state with high efficiency. Finally, a 
second fluorescence detection is applied for these atoms. 
In this way, we can determine the transition probability 
regardless of shot-to-shot atom number fluctuations. 

To lock the clock laser to the unperturbed ^Sq ~ "^Pq 
transition frequency, the Amp = 0, mp = ±9/2 lines are 
each interrogated on the high and low frequency points 
of maximum line slope i.e. at a detuning of about ±5 Hz 
from their respective line centres. The frequency shifts 
to address the Zeeman components are introduced by 
the AOM used to pulse the light of the clock laser. The 
magnetic field causes a splitting of the two lines of about 
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FIG. 3: Total fractional Allan deviation of the interleaved 
stabilization signal (dots, cumulated data set of all days), 
the frequency measurement against the H-maser (diamonds, 
one record), and of the cumulated measurement against CSFl 
(triangles). The line indicates the stability of the interleaved 
signal of 5 X 10"''^/ \/t/s, in which common mode effects 
do not show up. However, the uncertainty budget suggests 
that no effects occur down to the 10~^^ level that degrade 
the stability. The dashed line shows the stability of CSFl of 

1.4 X lo-^v^ATs- 

227 Hz. With the four excitation probabilities the clock 
laser can be steered to the mean of the two line cen- 
tres, which is in first order Zeeman-shift-free. The lock- 
ing is achieved by changing the frequency of the offset 
AOM between laser and cavity. From the observed tran- 
sition probabilities the Zeeman splitting is monitored and 
tracked, too. 

Under the conditions described above a single interro- 
gation cycle requires ~625 ms. We first interrogate one 
Zeeman component on both sides of its maximum and 
then switch to the second one. After these four interro- 
gation cycles we correct the determined frequency offset. 
To minimize locking errors, the drift of the clock laser due 
to the drift of its reference cavity is zeroed by a contin- 
uous frequency sweep on the offset AOM between laser 
and cavity. The frequency of the direct digital synthe- 
sizer (DDS) is altered in small steps (resolution 1.4 /xHz) 
at a rate of 100 Hz to obtain a slow quasi-linear ramp 
to compensate the initially estimated cavity drift rate. 
We update the sweep rate of the DDS after every inter- 
rogation of the four points. Thus the light of the clock 
laser is at a fixed frequency with a constant offset from 
the ^^Sr clock transition given by the switching AOM to 
the experiment. This light is delivered via an actively 
length stabilized optical fibre to a fs frequency comb, 
where the frequency of the clock laser is measured against 
a H-maser whose frequency is continuously measured by 
the PTE Cs fountain clock CSFl (ii,^. 

For systematic investigations we interleave two stabi- 
lization sequences with four interrogations in each case 
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30l 31 1 . With this method we measure systematic fre- 
quency changes due to different physical conditions in 
both sequences. This method has been apphed also dur- 
ing the actual frequency measurement with equal se- 
quences. The difference between the two stabilized fre- 
quencies provides information on possible locking errors 
and the quality of the stabilization of the clock laser 
on the atomic resonance. The observed fractional sta- 
bility of the difference frequency follows the dependence 
5 X / yjr/s (see Fig. which suggests a possible 

stability of 2.5 x 10~^^/ \/t/s if all interrogations are uti- 
lized for one laser stabilization. No significant frequency 
offset between the two interleaved stabilizations was de- 
tected. 
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SYSTEMATIC SHIFTS OF THE 
TRANSITION 



SR 



The interleaved stabilization method (see Sec. ^ allows 
to determine frequency shifts with a fractional statistical 
uncertainty of 10~^^ after 2 500 s (Fig. [J]). Similar uncer- 
tainty contributions have been evaluated in other groups 
with mostly comparable results d, [H, [l^, [2^ . The var- 
ious corrections and uncertainty contributions are listed 
in Tab. H] and will be discussed in detail below. 



ac Stark shifts 

Frequency shifts of the clock transition due to ac Stark 
shifts lead to the most prominent contributions to the 
corrections and uncertainties listed in Tab. HI The light 
shift of the lattice can be classified into first-order shifts 
with linear dependence on the trap depth, hyperpolar- 
izability contributions [3^ |33| (quadratic in trap depth) 
and contribution due to Ml or E2 transitions |33l. |34| . 

We have determined the (linear) Stark shift cancella- 
tion frequency of the lattice by variation of its frequency, 
using an interleaved stabilization with variation of the 
lattice depth between ^130 Er and '^350 Er- The ob- 
served frequency shifts are corrected for a small contri- 
bution from the hyperpolarizability [s^. The shifts are 
shown as function of the lattice laser frequency in Fig. 2] 
They have been normalized to the difference of the aver- 
age lattice trap depths experienced by the atoms. This 
averaging depends on the lattice geometry (deduced from 
the vibration frequencies) and atom temperature i.e. the 
mean particle energy in the potential [35|. The temper- 
ature is determined from the intensity ratio of red and 
blue axial sideband. The measured sideband frequen- 
cies were in agreement with the ones calculated using 
the measured power in the lattice and the polarizability 
at the magic wavelength 36, 33] . From this agreement 
we conclude that the lattice has a modulation depth as 
expected for perfect spatial overlap and that a contribu- 



FIG. 4: Clock transition shift measured with an interleaved 
stabilization between a trap depth of ~130 Er and ~350 Er 
(circles, upper panel). The line is a linear fit to the data to 
determine the zero crossing for \mF\ = 9/2. Lower panel: 
residuals of the fit with error bars. The zero crossing is found 
to be at a lattice laser frequency of 368 554 502(15) MHz. 



TABLE L Frequency corrections and their uncertainties 



Effect 



Correction Uncertainty 









lattice ac Stark (scalar & tensor) 


-0.2 


0.5 


hyperpolarizability 





0.2 


E2/M1 Stark shifts 





0.2 


ac Stark effect probe laser 


0.2 


0.2 


BBR ac Stark (room temperature) 


52.5 


1.3 


BBR ac Stark (oven) 


0.2 


0.2 


2nd order Zeeman 


0.28 


0.03 


collisions 


0.1 


0.2 


gravitation 


0.1 


0.1 


line pulling 





0.4 


servo error 





0.06 


AOM chirp 





0.2 


efliciency switch AOM 





0.02 


tunneling 





0.16 


2nd order Doppler 





< 10"^** 


total Sr 


53.2 


1.5 


fs comb and 200 m rf link 





1.5 


rf electronics 





4.0 


realization of the second 




7.6 


measurement statistics (60 500 s) 




5.7 


total all 


53.2 


10.5 



tion to the light shift from imperfect interference of the 
lattice beams can be neglected. 

It must be noted that the actual Stark shift cancella- 
tion frequency depends on the polarization of the light 
field, its orientation with respect to the magnetic field 
axis, and the Zeeman components considered [l^ l38| . 
We avoid possible vector light shift contributions due to 
the linear polarization of the lattice light field and its 
parallel orientation with the magnetic field axis plus the 
averaging of the transition frequencies of mp = ±9/2 
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components. We find the corresponding cancellation fre- 
quency to be at 368 554 502(15) MHz. The value is in 
agreement with the values given in 0, [3| • During the fre- 
quency measurement the lattice was operated at a 9 MHz 
lower frequency. We apply a correction of —2(5) x 10^^^ 
based on the mean lattice field seen by the atoms and the 
slope of the line in Fig. |4l We find a slope of the curve 
in Fig. 2] that is smaller than the slope determined by 
other groups [s^ . This inconsistency has however other- 
wise only a minor influence. The uncertainty reflects the 
inconsistency and the accuracy of the determination of 
the magic wavelength. 

The ac Stark shifts due to the hyperpolarizability and 
M1/E2 transitions are small [s^l. We use the trap depth 
of 125 Er during the frequency measurement to calculate 
an upper bound for these shifts. We deduce for both cases 
an uncertainty of 2 x 10^^^. 

The ac Stark effect of the probe laser is very small 
and can be calculated with sufficient accuracy. With a 
clock laser intensity of less than 0.6 mW/cm^ and the 
coeflacient from Ref. ^ of S^ock = 13(2) Hz/(W/cm2) 
we flnd an upper limit of 2(2) x 10"^^ for the correction. 
The beams of repump and cooling lasers are blocked by 
mechanical shutters and the vacuum chamber is shielded 
from remaining stray laser light. These measures allow 
neglecting possible ac Stark shifts by laser beams other 
than the ones from the clock and lattice laser. 

An important correction in the case of Sr frequency 
standards is due the blackbody radiation (BBR) shift. 
The shift close to room temperature is according to [s^ 



<5i^BBR -2.354(32) 



T 



300 K 



Hz. 



(1) 



To the uncertainty from theory, the uncertainty of the 
temperature measurement has to be added. We have 
continuously monitored the temperature of the vacuum 
system at ten positions with PtlOO sensors. Typical tem- 
perature variations over a day are at each point smaller 
than 100 mK, however peak-to-peak gradients of 4.2 K 
are observed over the vacuum system. The tolerance of 
the PtlOO sensors add an uncertainty of about 0.5 K. We 
calculate the temperature and its uncertainty assuming 
that it lies with equal probability an interval of 5.2 K 
around the best estimate of 296.8 K (rectangular prob- 
ability distribution, [4l|) and find a BBR correction of 
5.25(13) X 10-i5_ 

Further, we consider the influence of radiation coming 
from the hot furnace. Though it has no direct line-of- 
sight to the atoms, blackbody radiation could be multiply 
scattered in the water-cooled Zeeman slower tube, enter 
the vacuum system and be reflected onto the atoms. In a 
worst case scenario this can be approximated by a virtual 
orifice placed on the wall of the vacuum system that has 
the area of the collimation pinhole shielding the furnace 
BBR (Sec.p) and that emits BBR of 793 K homogeneously 
into the vacuum system. 



Using the models described in Ref. [23'| this radiation 
causes additional shifts by the direct line-of-sight between 
orifices and atoms of 



'54br(793 K) « X 146 Hz 



and by the diffusely scattered radiation of 



(2) 



<^4br (793 K) « X 146 Hz. (3) 

47r a 



e 1 



We have introduced 9 as the solid angle under which 
the atoms see the virtual orifice and a, the absorption 
coefficient of the walls. Using the diameter of our vac- 
uum system of 335 mm and assuming a « 0.2 [4^ . we 



find 



;dir) 
BBR 



< 1 mHz and 



(diff) 

BBR 



< 7 mHz. These 



corrections are on the level of 2 x 10 



Zeeman effect 

As the laser is locked to the average frequency of the 
Amp — Zeeman transitions from levels mp — -1-9/2 
and mp — —9/2 the linear Zeeman shift is removed for 
constant fields. Drift in the field give rise to a residual 
first order Zeeman shift. Great care has been employed 
that the magnetic field is the same when both Zeeman 
components are probed. Nevertheless small drifts of the 
offset field have been observed as temporal variations of 
the frequency splitting of the two interrogated Zeeman 
lines. They cause a frequency drift of the two Zeeman 
lines in opposite direction, which is less than 0.23 mHz 
per single loading and interrogation cycle. This causes a 
frequency offset of only 0.23 mHz after the four interro- 
gations because the effect is mostly compensated by the 
opposite frequency change. We use this value as uncer- 
tainty for the residual linear Zeeman shift and have an 
uncertainty of 5 x 10""'^^, which we neglect. 

The stabilization sequence provides a value of the Zee- 
man splitting, which is needed to correct for the sec- 
ond order Zeeman shift. We have varied the homoge- 
neous magnetic field during the spectroscopy pulse be- 
tween 23 /iT and 630 /zT and have measured the Zeeman 
shift with the interleaved stabilization against a refer- 
ence of 23 fxT. The magnetic field was calibrated using 
the linear Zeeman shift coefficient given in . Figure [5] 
shows the measured shifts in the upper panel together 



with a fit of the type aB 



In the lower panel the 



fit residuals are shown. We find for the quadratic coeffi- 
cient a — —23.0(3) Hz/mT^ which is in good agreement 
with the value from Ref. [38,]. For the magnetic field 
of 23(1) /^T during the frequency measurements we ob- 
tain a correction for the second order Zeeman shift of 
0.28(3) X 10-16. 
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FIG. 5: Measurement of the second order Zeeman shift (top, 
dots) with a parabola fit of the type aB^ + c and residuals 
(bottom). 



Density shift 



-0.2 - 



1000 



2000 



3000 



4000 



5000 



population difference (arb. units) 



FIG. 6: Measurement of the density shift (dots) and linear fit 
without an offset (line) and its statistical uncertainty (dashed 
lines). The frequency measurements were performed at a 
population of about 600 in these arbitrary units, which cor- 
responds to about 10* atoms. The error bars indicate the 
statistical uncertainty of the interleaved stabilization signal 
of the corresponding measurement. 



Observations of a density shift with spin polarized 
fermions and their theoretical interpretation have trig- 
gered intense discussions SSSH. A collisio n shift 
is observed if different atoms experience different Rabi 
excitations due to e.g. spatial inhomogeneous excitation 
fields or motional state dependent Rabi frequencies. Care 
has been taken to avoid these issues in our experiment 
by careful alignment and high beam quality. To investi- 
gate a possible density shift we vary the atom number in 
the lattice by a variation of the length of the first MOT 
phase. We have not observed a significant dependence 
of the temperature in the lattice on the atom number 
and therefore assume the density to be proportional to 
the atom number. In our setup, the determination of 
absolute atom numbers can only be done with large un- 
certainties with ®''Sr, but relative atom numbers can be 
determined reproducibly with high precision. We have 
increased the atom number by nearly a factor of ten 
compared to the conditions during the frequency mea- 
surements and have measured a possible collision shift 
with respect to a reference density comparable to the 
one of the frequency measurement. The data are plotted 
in Fig. [6l Theory suggests [13] that a density shift may 
sensitively depend on the excitation parameters which 
contributes to the scatter of the data. The uncertainty 
of the fitted slope shown in Fig.lHlreflects this effect. Dur- 
ing the frequency measurement the excitation probability 
was not as well controlled as during this measurement of 
the coUisional shift. Thus, for the correction of the fre- 
quency measurement we have increased the uncertainty 
of the fitted slope by a factor of three and deduce a cor- 
rection due to collisions of 0.1(2) x 10""'^^. 



Other corrections and uncertainties 

Further corrections have to be applied: The effective 
position of the atoms in the Cs fountain is 11(10) cm 
above the position of the strontium atoms, which leads 
to a correction of 0.1(1) x 10^^^ due to the gravitational 
red shift. The height difference was determined by tri- 
angulation between the two buildings that are separated 
by about 200 m. 

The magnetic field applied during spectroscopy causes 
a Zeeman splitting of ^25 Hz between neighbouring lines. 
Line pulling can appear due to atoms in undcsired Zee- 
man levels that are not removed during the preparation 
sequence. From the achieved signal quality (Fig. [2]) we 
can rule out a population fraction in these levels of larger 
than 5%, which leads to an uncertainty due to line pulling 
of 4 X 10-1^. 

The steering of the clock laser to the atomic resonance 
frequency is compromised by a cavity drift that varies in 
time. Locking errors occur if the cavity drift rate changes 
too rapidly. We have observed changes of the cavity drift 
rate Dcav of up to 2 x 10^^ Hz/s^. The lock of the laser 
to the reference line employs two parts. First, a fraction 
G = 0.7 of the measured frequency offset z/err is applied 
directly as a correction. Second, the drift rate is altered 
by /^I'orr/Tcycio with /3 = 0.05. The time constant for the 
settling of the drift rate correction is 



Tdrift 



'^cyclc 



G 



ln(l-| 



'cycle ~ 



(4) 



which amounts to 35 s for our cycle time Tcycio = 2.5 s. 
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For a cavity with a constant change of the drift rate of 
i'cav the frequency error z^iockcn- is determined based on 
the assumption that in steady state the correction of the 
drift rate has to match the change in the cavity drift that 
occurs within that lock cycle: 



cav ■ '^cyclc 



"^cycle 



^lock err 



^lock err ^c. 



'^cyclc 



(5) 

For the given experimental conditions this locking error 
leads to a very small uncertainty of 6 x 10~^*. 

We mentioned in Sec. □ that we have implemented a 
fibre length stabilization that operates in a pulsed mode. 
It is used to correct for phase errors due to expansion 
of the fibre, AOM effects [s^l, movement of the lattice 
sites holding the atoms with respect to the clock laser, 
or other effects that change the interrogation phase of 
the clock laser at the position of the atoms. However, 
the servo loop can be activated only at the beginning of 
the spectroscopy pulse when an error signal is available. 
This leads to a phase transient before the servo loop set- 
tles that causes phase chirps at the position of the atoms. 
We have measured these chirps with a heterodyne detec- 
tion and have calculated the resulting frequency offsets 
using a numerical integration of the time evolution of 
the density matrix [27]. The remaining uncertainty due 
to phase errors is 2 x 10^^^. 

The frequency detuning of the clock laser to address 
the Zeeman components is accomplished by the switching 
AOM frequency by about 5 Hz. This may result in a sys- 
tematic variation of the laser intensity and thus systemat- 
ically different excitation probabilities at the lock points. 
By measuring the light power transmitted through the 
fibre delivering the clock laser light to the atoms as a 
function of AOM frequency we found that this introduces 
only a marginal uncertainty of 2 x 10^^*. 

Many Sr experiments use vertically oriented optical 
lattices to avoid tunneling in shallow traps by shifting 
of the lattice sites by the gravitational potential [45| . 
We have chosen a horizontal orientation to achieve bet- 
ter spatial overlap between MOT and lattice and thus 
to load a larger number of atoms into the lattice. This 
requires however a deeper potential to avoid tunneling 
between the sites and the connected broadening of the 
band structure. 

With the experimental parameters 77% of the atoms 
populate the two lowest bands of the lattice. We use the 
half width of the first excited band of the lattice as upper 
bound for the effect of tunneling leading to an uncertainty 
contribution of 1.6 x 10"^^. 

As last uncertainty contribution we consider the second 
order Doppler effect of the atoms in the lattice. This 
effect is however <Si 10~^^ and completely negligible. 

All corrections discussed in this section add up to 
53.2 X 10^^^ with an uncertainty of 1.5 x 10"^^ (see 
Tab.m. 



FREQUENCY MEASUREMENT AND RESULTS 

The frequency of the PTE *^Sr frequency standard was 
measured on three days in October 2010. The frequency 
of the clock laser was measured with a 1.5 /im Er-doped 
fs fibre frequency comb referenced to a 100 MHz rf sig- 
nal delivered by a semi-rigid cable with low temperature 
sensitivity from a H-maser. The 5 MHz output of the 
H-maser was continuously measured by the Cs fountain 
clock CSFl. The stability (indicated by the total Allan 
deviation [46j) between the Sr frequency standard and 
the H-maser or fountain clock is shown in Fig. |3l 

Both beat frequencies from the Sr clock laser with the 
fs comb and of the / — 2/ interferometer generating the 
carrier envelope offset frequency were tracked twice by 
tracking oscillators. The signals of all tracking oscilla- 
tors were counted as was the repetition rate of 100 MHz. 
The double tracking was applied to identify rare cycle 
slips of phase locked loops and discard them from the 
data set. Frequency comb measurements connecting mi- 
crowave and optical frequencies have been demonstrated 
by several groups at an accuracy level of 10~^^ or be- 



low [471 l48l for averaging times exceeding 1000 s. In the 
measurements shown here, however, the rf-reference was 
delivered to the frequency comb via a 200 m semi-rigid 
cable, which connects the Cs-fountain clock system (in 
building A) with the frequency comb and Sr-frequency 
standard (in building B). We thus investigated this 200 m 
cable with respect to the phase stability of the deliv- 
ered 100 MHz reference signal; at an averaging time of 
10 000 s, we observe frequency fluctuations (Allan devi- 
ation) of about 1.5 X 10~^^. 

To give an upper bound for the uncertainty arising 
from the frequency comb measurement, including the 
cable connection, we additionally performed simultane- 
ous measurements of an optical frequency (from a cav- 
ity locked laser in building B) versus an rf-reference (a 
H-Maser in building A) with two frequency combs lo- 
cated in the two different buildings. Measurements by 
the two fs-combs agreed at the level of statistical uncer- 
tainty of 1.5 X 10^^® for an averaging time of 10 000 s, 
which suggests that they were limited in stability by the 
rf-link discussed above. During the three days when the 
Sr-clock transition frequency was measured the air con- 
ditioning system in the room of the fountain clock gave 
rise to periodic 0.2 K peak-to-peak temperature varia- 
tions with periods in the range of 10 000 s. An investi- 
gation of thereby induced phase shifts of the fountain rf 
electronics results in an uncertainty contribution of less 
than 4 x 10^^® taking into account the periodicity of the 
variations and the durations of the single measurements. 

The Cs fountain contributes with a systematical uncer- 
tainty of 7.6 X 10^^^. The fountains instability of 1.4 x 
10""'"^/ \/t/s and the 60 500 s total record length of the 
three days measurement lead to a statistical uncertainty 
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FIG. 7: Frequencies of the *''Sr clock transition measured by 
different laboratories: Paris (square, [1^), Boulder (circle, 
[19I]'). Tokyo (triangle, [l3|)i and Braunschweig (diamond). 
The vertical line gives the recommendation for *^Sr as sec- 
ondary representation of the second jl^] with its uncertainty 
(dashed lines). 



of 5.7 X 10 . The measured frequencies of all days are 
in agreement with each other. The weighted average of 
the three frequency values is 429 228 004 229 872.9(5) Hz. 
It is i n g ood agreement with previous measurements 
0, [13, Il9l. I20I and the frequency value recommended for 



the secondary representation of the second by ^''Sr 
(see Fig.IZl). 



13| 



the magic wavelength during clock operation than was 
possible during the present measurement could reduce 
the biggest contribution to the lattice ac Stark effect in 
Tab. H] below 0.1 x 10~^®. The uncertainty from the line 
pulling effect due to atoms in undesired Zeeman levels 
(mainly the mp — ±7/2 levels) can be reduced by a 
more detailed investigation of its influence in the actual 
preparation conditions by employing different magnetic 
fields with the associated different influence on the line 
pulling effect. Similarly, there is no principal limitation 
that should prevent to reduce the other contributions to 
below 10~^^ thereby leading to a combined fractional un- 
certainty around 10~^^. With several ^^Sr clocks being 
investigated in different institutes it is expected that an 
in-depth evaluation of this type of optical clock on the 
lO"^*" level will be achieved soon. Having different opti- 
cal clocks based on neutral atoms and on single ions at 
hand that have an order of magnitude higher accuracy 
than the best primary Cs clocks of today, the next steps 
towards a re-definition of the second as the base unit of 
the SI could be made. 
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CONCLUSION 



We report on a frequency measurement of the ^Sq - 
■^Fq clock transition in ®^Sr with a fractional uncertainty 
of 1 X 10""'^^ limited by the comparison with the Cs atomic 
clock. The measurement is in good agreement with previ- 
ous results achieved in three other laboratories within the 
combined uncertainties. The estimated standard uncer- 
tainty of our ^^Sr standard by itself was estimated to be 
below 2 X 10^^® which competes well with the uncertainty 
of the best primary Cs clocks for the realization of the 
SI units s and Hz and, hence, supports the choice of this 
standard as secondary representation of the second. The 
reported uncertainty represents the status of the clock 
during the last frequency measurement and could be fur- 
ther improved. 

The presently largest contributions to the uncertainty, 
results from the blackbody ac Stark effect, the lattice ac 
Stark effect, and the line pulling effect with 1.3 x 10^^^, 
0.5 x 10^^^, and 0.4 x 10~^^, respectively. Investiga- 
tions are under way to determine the influence of the 
blackbody radiation shift 23| with the prospect to re- 
duce the associated correction and uncertainty signifi- 
cantly. Tuning the frequency of the lattice laser closer to 
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